Quantum phase transitions in a two-dimensional quantum XYX model: 
Ground-state fidelity and entanglement 
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A systematic analysis is performed for quantum phase transitions in a two-dimensional anisotropic 
spin 1/2 anti-ferromagnetic XYX model in an external magnetic field. With the help of an inno- 
vative tensor network algorithm, we compute the fidelity per lattice site to demonstrate that the 
field-induced quantum phase transition is unambiguously characterized by a pinch point on the 
fidelity surface, marking a continuous phase transition. We also compute an entanglement estima- 
tor, defined as a ratio between the one-tangle and the sum of squared concurrences, to identify 
both the factorizing field and the critical point, resulting in a quantitative agreement with quantum 
Monte Carlo simulation. In addition, the local order parameter is "derived" from the tensor network 
representation of the system's ground state wave functions. 
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Quantum critical phenomena are crucial in our under- 
standing of the underlying physics in quantum many- 
body systems, especially in condensed matter systems, 
due to their relevance to high- Tc superconductors, frac- 
tional quantum Hall liquids, and quantum magnets 
The latest advances in this area arise from quantum infor- 
mation science. Indeed, various entanglement measures 
have been widely applied to study condensed matter sys- 
tems. Remarkably, for one-dimensional (ID) quantum 
systems, von-Neumann entropy, as a bipartite entangle- 
ment measure, turns out to be a good criterion to judge 
whether or not a system is at criticality 0, 0, @, d, 0j liQ- 
On the other hand, fidelity, another basic notion in quan- 
tum information science, has demonstrated to be funda- 
mental in characteri zing phase transitions in quantum 
many-body systems [lOl . [ill . Il^ . This adds a new 
routine to explore quantum criticality in condensed mat- 
ter physics from a quantum information perspective. 

However, with only a few notable exceptions fl3. 
not much work has been done for two-dimensional (2D) 
quantum systems, due to great computational challenges. 
In fact, despite the existence of well-established numer- 
ical algorithms, such as exact diagonalization, quantum 
Monte Carlo (QMC), the density matrix renormalization 
group (DMRG) and series expansions, drawbacks become 
obvious when one deals with frustrated spin systems. A 
typical example is the QMC, which suffers from the no- 
torious sign problem. However, a promising progress, 
inspired by new concepts from quantum information sci- 
ence, has been made in classical simulations of quantum 
many-body systems. The algorithms are based on an 
efficient representation of the system's wave functions 
through a tensor network. In particular, matrix prod- 
uct states (MPS) 1^ 17, [l3|, a tensor network already 



jected entangled-pair states (PEPS) constitute the basis 
to simulate 2D quantum lattice systems [2ll . [2^ . 

The aim of this paper is to show that the fidelity per 
lattice site, first introduced in Ref. [ll|, is able to un- 



veil quantum criticality for a 2D anisotropic spin 1/2 
anti-ferromagnetic XYX model in an external magnetic 
field. This is achieved by exploiting tensor network algo- 
rithms, i.e., innovative algorithms inspired by the latest 
achievements in our understanding of quantum entan- 



glement [21|, [2^ . We show that the field-induced quan- 
tum phase transition is unambiguously characterized by a 
pinch point on the fidelity surface, marking a continuous 
phase transition. In addition, we compute an entangle- 
ment estimator, defined as a ratio between the one-tangle 
and the sum of squared concurrences (for all the pairwise 
entanglement between two spins), to identify both the 
factorizing field and the critical point, resulting in con- 
sistent conclusions as drawn from the fidelity approach, 
with an extra result about a factorizing field hf. Our 
results are compared to those of the QMC simulation 
by Roscilde et al. for the model. We stress that 
the QMC simulation is carried out for a system on a fi- 
nite square lattice at very low but finite temperatures, 
whereas our simulation is directly performed for an infi- 
nite system at zero temperature. 

Quantum XYX model. We consider the 2D antiferr- 
magnetic spin-1/2 XYX model in a uniform z axis exter- 
nal magnetic field: 



H = J 



<i,j> 



( QX nx 



StS^^) + Y,hSl (1) 



present in DMRG, are used in the time-evolving block 
decimation (TEBD) algorithm to simulate time evolu- 
tion in ID quantum lattice systems 1^, 23], whereas pro- 



where J > is the exchange coupling, < i,j > runs 
over all the possible pairs of the nearest neighbors on a 
square lattice, and h is the external magnetic field. From 
the non-commutativity of the spin-1/2 Pauli operators, 
XYX model is expected to undergo a continuous quan- 
tum phase transition, with the same universality class as 
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the 2D quantum Ising model in a transverse field. < 1 
and Aj, > 1 correspond to easy-plane (EP) and easy-axis 
(EA) behaviors, respectively. The ordered phase in the 
EP (EA) case arises from spontaneous symmetry break- 
ing along the x (y) direction, with a finite value of the 
order parameter, i.e., the magnetization rrix (jny) below 
the critical field he- In Ref. [2^, the QMC simulation 
was exploited to discuss the connection between quantum 
phase transitions and entanglement measures, where an 
entanglement estimator, defined as the ratio between the 
one-tangle and the sum of squared concurrences, was sys- 
tematically analyzed to signal a quantum critical point. 

The infinite projected entangled-pair state (iPEPS) 
method. Let us briefly recall the iPEPS algorithm. Con- 
sider a finite two-dimensional square lattice where each 
site, labeled by a vector r ~ ix,y), is represented by a 
local Hilbert space V^^ = C"^ of finite dimension d. Let 
a vector j^*) denote a pure state in the (global) Hilbert 
space and the operator H = X^rr' h^^^^ be a Hamilto- 
nian with the nearest neighbor interactions on the lattice. 
Each lattice site has been represented by a tensor A^^, so 
a PEPS for the state l^*) consists of a set of tensors A^^ . 
The tensor A^^^^^^ is made of complex numbers labeled 
by one physical index s and four inner indices u, d, I and 
r. The physical index runs over a basis of V^^, so that 
s = 1, • ■ • , d, whereas each inner index takes D values, 
where D is some bond dimension, and connects the ten- 
sor with the tensors in the nearest neighbor sites. Thus, 
in a lattice with N sites, a PEPS depends on 0{ND^d) 



parameters 21 1. 



Now we move to a system defined on an infinite square 
lattice and assume that both |^) and H are invariant 
under shifts by two lattice sites. We exploit this in- 
variance to store the iPEPS using only two different 
tensors A and B. Given an iPEPS for a state 
(e.g., a product state), the iPEPS algorithm allows to 
perform an evolution in imaginary time to compute a 
ground state wave function of a g iven Hamiltonian H , 
|*,) = e-«-|vI/o)/||e-^-|x],o)|| [if. 

A contraction process, which is related to an evolution 
task, is done in order to get the effective environment 
for a pair of tensors A and B 



22l |. In practice, a global 



optimization problem has been reduced to a local two 
site optimization problem. Two new tensors A' and B' 



could be computed by a sweep technique [17|, originally 
devised for an MPS algorithm applied to ID quantum 
systems with periodic boundary conditions [l8l |. 

The fidelity per lattice site. Consider a finite 2D square 
lattice system described by Eq. dT]), with the external 
magnetic field /i as a control parameter. For two differ- 
ent ground states of the system, '^{h) and ^{h'), corre- 
sponding to two different values h and h' of the control 
parameter, respectively, the fidelity per lattice site d is 
defined as: 

\nF{h,h') 




FIG. 1: (color online) The fidelity per lattice site d{h,h'), 
as a function of h and h' for two ground states of the two- 
dimensional quantum XYX model. This defines a 2D fidelity 
surface embedded in a three-dimensional Euclidean space. A 
continuous phase transition point he « 3.489 is characterized 
as a pin ch point {hc,hc) on the fidelity surface, as argued in 
Ref. Here we have taken the bond dimension D = 2. 

The gray line denotes the normalization: d{h,h) = 1. 



where N is the system size, and F{h, h') = | {4'{h')\'ip{K)) \ 
is the ground state fidelity. The fidelity per lattice site 
d depicts how fast the fidelity goes to zero when N gets 
large. Remarkably, the fidelity per lattice site d is well 
defined in the thermodynamic limit: 



In d(h, h') = lim 



\iiF{h,h') 
N ■ 



(3) 



\Yid{h,h') 



N 



(2) 



It satisfies the properties inherited from fidelity F{h, h'): 
(i) normalization d{h,h) = 1; (ii) symmetry d{h,h') = 
d[h' , h)- and (in) range < d{h, h') < 1. 

As shown in Ref. [ll[, the fidelity per lattice site 
d{h,h') succeeds in capturing nontrivial information 
about stable and unstable fixed points along renormaliza- 
tion group flows. Specifically, the fidelity surface, defined 
by the fidelity per lattice site d{h, h') as a 2D surface em- 
bedded in a three-dimensional Euclidean space, exhibits 
singularities when h — he or h' — he. That is, d{h,h') 
exhibits singular behaviors when h crosses he for a fixed 
h', or h' crosses he for a fixed h. Therefore, a phase 
transition point he is characterized as a pinch point (2^ 
(he, he) for continuous QPTs, i.e., the intersection of two 
singular lines h = he and h' — he. 

The fidelity per lattice site may be computed from 
the iPEPS representation of the ground state wave func- 
tions, following the transfer matrix approach described 
in Ref. [14|. We plot d{h,h') in Fig. [TJ computed with 
the help of the iPEPS algorithm [2^ with bond dimen- 
sion D = 2 (the result for 13 = 3 is very similar to that 



3 



for D = 2). A pinch point on the fidehty surface defined 
by d{h, h') as a function of h and h' clearly indicates a 
second order phase transition. In addition, the two sta- 
ble fixed points at ft. = and h = oo are characterized as 
the global minima of the fidelity surface (for a fixed Aj^). 

The one-tangle and the concurrence. We now exploit 
the iPEPS algorithm to extract the ground state entan- 
glement properties of the quantum XYX model on an 
infinite square lattice. We first compute the one-tangle 
Ti, defined as Ti = 4detp^^\ where p'-^' is the single site 
reduced density matrix, for a specific value Aj^ = 0.25 
with D — 2 and 3. However, we stress that, our dis- 
cussion, although only confined to the Ay = 0.25 case, 
is actually quite generic and applies to all other values 
of Aj, < I. The one-tangle reflects the entanglement 
between a single site and the rest of the system. Note 
that there exists a factorizing field, at which the one- 
tangle Ti vanishes j2^. The exact theoretical value of 
the factorizing field /i/ = 2^2(1 + A^) for Ay = 0.25 
is approximately 3.162, consistent with the iPEPS re- 
sults up to four digits. When the external field is in- 
creased beyond the factorizing field ft/, a cusp occurs for 
an entanglement ratio i? = T2/T1 at a critical point ftc, 
where T2 denotes the sum of squared concurrences (for 
all the pairwise entanglement between two spins). The 
cusp in the entanglement ratio can be regarded as a sig- 
nal of a quantum phase transition. In Fig. [21 both the 
one-tangle ti and the sum of squared concurrences T2, 
obtained from the iPEPS with D = 3, are greater than 
those for D = 2. This is due to the fact that for larger D, 
the PEPS representation accommodates more entangle- 
ment. The critical point ftc ~ 3.485 from the iPEPS for 
L» = 3 is smaller than that for £> 2 (ft^ « 3.489). The 
shift is quite small, so one may expect that the critical 
point for small D does not deviate significantly from the 
converged result with large D. Our result indicates that 
the iPEPS algorithm is able to capture the entanglement 
properties, and gives rise to the results consistent with 
the QMC simulation. 

Local order parameter. The efficient tensor network 
representation of the system's ground state wave func- 
tions makes it possible to extract an (optimized) local 
order parameter, according to a general scheme advo- 
cated in Ref. [l^- In fact, once the critical field ftc is 
determined, one may choose two representative ground 
states, one for an external magnetic field ft less than the 
critical field ftc and the other for an external magnetic 
field ft greater than the critical field ftc. Then the re- 
duced density matrix p^^^ for a single lattice site in an 
infinite-size lattice is computed for two different values 
of the external magnetic field ft, corresponding to ft > ftc 
and ft < ftc, respectively. It is readily found that the 
one-site reduced density matrix p^^-* displays different 
nonzero-entries structures in two phases, with (Sx) be- 
ing zero for h > he and nonzero for h < he. Also note 
that the Z2 symmetry is spontaneously broken, since the 
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FIG. 2; (color online) Upper panel: The one-tangle ri and 
the sum of squared concurrences T2 as a function of an ap- 
plied external field h for the 2D quantum XYX model with 
Ay — 0.25. The data are presented for both D — 2 and 
D = 3. The factorizing field, at which the one-tangle ri van- 
ishes, is indicated by an arrow labeled by ft/ around 3.162. 
Above the factorizing field ft/, a steep increase of ri and T2 re- 
flects a rapid increase of entanglement around a critical point 
he ~ 3.489 {D = 2) and he ~ 3.485 (D = 3). Lower panel: 
An entanglement ratio R — T2/T1, exhibits a cusp, a signal of 
a continuous phase transition, at a critical point he obtained 
using the infinite projected entangled-pair state (iPEPS) al- 
gorithm with D = 2 and 3. This is in an agreement with that 
resulted from the fidelity approach. 

reduced density matrix p^^^ does not commute with the 
symmetry generating operator in the symmetry-broken 
phase h < he- This implies the existence of a local or- 
der parameter: rUx = {'4>(h)\Sx\'4>{h)), characterizing the 
second order phase transition that belongs to the same 
universality class as that of the 2D quantum Ising model 
in a transverse field. 

Another interesting local observable rriz = 
{il>{h)\Sz\ip{h)) is also carefully investigated, which 
exhibits singularities at the critical field ftc. If the 
external transverse magnetic field ft is raised from 
below the critical field ftc, then the z magnetization is 
a monotonic curve and gradually reaches a saturated 
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locate it, as long as such a factorizing field exists for a 
system considered [27 1. 

This work is supported in part by the National Natural 
Science Foundation of China (Grant Nos. 10774197 and 
10874252). 



FIG. 3: (color online) The x magnetization mx{h) in the 
ground state \'^{h)) of the 2D quantum XYX model {Ay — 
0.25), a local order parameter, is readily read off from the in- 
finite projected entangled-pair state (iPEPS) representation 
of the ground state wave functions with D = 2 and 3. Dashed 
lines are a guide to the eye. For the external magnetic field h 
less than the critical magnetic field he, the order parameter 
takes a non-zero value and decays as h increases. Once the 
field is greater than the critical field he, the order parame- 
ter becomes zero. Notice that, around he, the derivative of 
the z magnetization mz{h) changes rapidly. Thus it exhibits 
singularities at the critical field he- 



value. 

We plot rux and ruz as a function of the external 
magnetic field h for the 2D quantum XYX model with 
Ay = 0.25 in Fig. [31 We have presented the data for both 
D = 2 and 3, with the truncation dimension x for IMPS 
used in the contraction of the iPEPS representation up 
to 30. On the other hand, although the factorizing field 
can be readily located in Fig. [21 there are no unusual 
features appearing in rrix and rUz around this point. 

Conclusions. We have performed a systematic analy- 
sis of both the ground state fidelity and entanglement for 
the 2D quantum XYX model. This is achieved by com- 
puting the ground state wave functions by means of the 
iPEPS algorithm. The results are compared with those 
obtained from the QMC simulation [23[, giving rise to 
a remarkable quantitative agreement for the factorizing 
field. From a fidelity perspective, a connection between a 
pinch point on the fidelity surface and a continuous phase 
transition point has been demonstrated, thus allowing us 
to determine the ground state phase diagram of the 2D 
quantum XYX model. Once this has been done, one 
may read off the local order parameter from representa- 
tive states in terms of the iPEPS representation of the 
ground state wave functions. 

An interesting point worth to be mentioned is that, 
in contrast to entanglement measures, the fidelity per 
lattice site fails to locate a factorizing field. However, 
one may resort to a closely related quantity, i.e., the so- 
called geometric entanglement introduced in Ref . [1^ , to 
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